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ABSTRACT 



Context. Members of nearby moving groups are promising candidates for the detection of stellar or substellar components by direct 
imaging. Mass estimates and magnitude measurements of detected companions to young stars are valuable input data to facilitate the 
refinement of existing pre-main-sequence stellar models. In this paper we report on our detection of a close companion to HD 160934, 
a young active star, SB 1 spectroscopic binary and a suggested member of the AB Doradus moving group. 

Aims. We obtained high angular resolution images of nearby young stars, searching for close companions. In the case of HD 160934, 
direct imaging is combined with unresolved photometry to derive mass estimates. 

Methods. High angular resolution was achieved by means of the so-called "Lucky Imaging" technique, allowing direct imaging 
close to the diffraction limit in the SDSS z' band. Our results are combined with pre-discovery HST archive data, own UBV(RI)c 
broadband photometry, published JHK magnitudes and available radial velocity measurements to constrain the physical properties of 
the HD 160934 close binary. 

Results. At an assumed age of ~80Myr, we derive mass estimates of 0.69 Mq and 0.57 Mq, respectively, for HD 160934 and its 
close companion. We suggest that the direct detection may be identical to the spectroscopically discovered companion, leading to a 
period estimate of ~8.5 years and a semimajor axis of a«4.5 A.U. 
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1. Introduction 

HD 160934 (= HIP 86346) is a y oung l ate-type star with a spec- 
tral type of K7 to M0 (Reid et al. 1995; Zuckerman et al. l2004al l 
at a distance of *24.5pc (Perryman et al. I19971 I. It is chromo- 
spherically active (Mulliss & Bopp |199"4l ) with prominent EUV 
(e.g., Pounds et al. 1993) as well as X-ray emission with an X- 
ray luminosity of L x = 3.4x 10 22 W (Hiinsch et al. [19991 The ac- 
tivity can also be traced in the Ha line, which is seen in emission 
with an equivalent width between -0.09 and -0.13nm (Mulliss 
& Bopp [T9941 Gizis et al. 120021 Zuckerman et al. I2004ab . The 
detection of the Li 6708 A line with an equivalent width of 40 mA 
(Zuckerman et al. 2004a ) gives further evidence that HD 160934 
is a relatively young star. These youth indicators combined with 
the 3d space motion led Zuckerman et al. (2004a 2004b} and 
Lopez-Santiago et al. (120061 ) to suggest that HD 160934 might 
be a member of the *50Myr old AB Dor moving group. 

Because of its proximity to the Sun and its young age, 
HD 160934 is a good candidate for the direct detection of 
substellar, or even planetary mass companions. McCarthy & 
Zuckerman ( 2004 ) report that no brown dwarf companion could 
be found at projected separations larger than 75 A.U. as a re- 
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Hispano Aleman (CAHA) at Calar Alto, operated jointly by the 
Max-Planck-Institut fiir Astronomie and the Instituto de Astrofisica 
de Andalucfa (CSIC). Based on observations obtained with MPIA's 
Konigstuhl mountain 70cm telescope. 



suit of a near infrared coronagraphic study carried out at the 
Keck observatory. Using HST/NICMOS in coronagraphic mode, 
Lowrance et al. (2005 ) report the detection of a possible wide 
companion to HD 160934 at a projected separation of w8'.'7 (cor- 
responding to ^210 A.U.) and at a position angle of «235 . The 
brightness difference between the companion candidate, desig- 
nated HD 160934B, and HD 160934 A is AH = 9.2 mag. Under 
the assumption that HD 160934 B constitutes a physical com- 
panion to HD 160934, Lowrance et al. (2005) derive a mass es- 
timate of »0.15 Mq for this companion. 

By combining 37 radial velocity measurements, Galvez et 
al. ([2006 ) were able to identify HD 160934 as a spectroscopic 
SB 1 binary and suggested a period of P=6246.23 1 8 days, a high 
eccentricity of e=0.8028, and a spectral type of M2-M3V for the 
close companion, so that HD 160934 may be actually at least a 
triple system. 

In July 2006, we started a high angular resolution survey for 
close stellar and substellar companions to young nearby stars, 
taking advantage of the diffraction limited performance facil- 
itated by using 'Lucky Imaging' at the Calar Alto 2.2 m tele- 
scope. The absence of a coronagraphic mask in our set-up, and 
the fact that the diffraction limit of a 2.2 m telescope at the effec- 
tive observation wavelength of 900 nm corresponds to an angular 
resolution of wC'l means that our survey is sensitive to com- 
panions at close separations. In this Research Note we report 
on the first results of the survey, namely the direct imaging of 
HD 160934 and a close companion. We summarise the available 
photometric measurements, give new UBV(RI)cz' magnitudes, 
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Table 1. Direct imaging observing log for HD 160934. 



Date 


Telesc./Inst. 


Filter 


tint 


June, 30 1998 


HST/NIC2 


F165M 


2 x 0.626s 


July, 8 2006 


CA 2.2m/AstraLux 


RG780 


4.5s° 


July, 8 2006 


CA 2.2m/AstraLux 


RG830 


4.5s" 



a 2% best of 15000 frames with an individual t CXD = 0.015s 



and discuss the possibility that the spectroscopic and our direct 
detection refer to the same companion. 

2. Observations and data reduction 

2.1. Direct imaging with AstraLux at the Calar Alto 2.2m 
telescope 

Observations at the Calar Alto 2.2m telescope were obtained 
with the new 'Lucky Imaging' instrument AstraLux. The Lucky 
Imaging technique is based on the fact that in a large set of 
ground based short exposure images at least a small fraction will 
be relatively unaffected by atmospheric turbulence. By finding 
and using only these high-quality images, one is able to obtain 
diffraction limited images from the ground without the need of 
adaptive optics (Baldwin et al. l2001l Tubbs et al. I20021 Law et 
al. l2006l 

The AstraLux instrument basically consists of an electron 
multiplying, thinned and back-illuminated CCD, which is oper- 
ated in frame transfer mode and capable of running at a frame 
rate of 34 Hz when reading the full 512x512 pixel frame. The 
camera is mounted at the Cassegrain focus of the telescope 
behind a simple Barlow lens assembly, which changes the fo- 
cal ratio from f/8 to f/32 in order to get an image scale of 
46.6 mas/pixel at a physical pixel size of 16xl6jt/m 2 . This pro- 
vides Nyquist sampling at A w lyum, and slightly undersamples 
the theoretical point spread function (PSF) at the typical obser- 
vation wavelengths of 700-950 nm. With the exception of the 
mount needed to attach the camera to the telescope, all compo- 
nents were bought off-the-shelf, making AstraLux a simple and 
cost-effective instrument, taking less than six months from the 
first idea to first light. 

HD 160934 was observed as part of a larger sample of nearby 
young stars on July 8, 2006 under good atmospheric conditions 
with a photometric sky and a V-band seeing around 078. The 
raw data for HD 160934 consists of two FITS cubes of 15000 
frames with 15ms exposure time each. We used only a 256x256 
pixel subarray of the chip, leading to a frame rate of ^65 Hz. The 
target was observed through two different filters, namely Schott 
RG780 and RG830 longpass filters, the latter roughly matching 
the SDSS z'-band (for the definition of the SDSS passbands, see 
e.g. Fukugita et al. 1996). 

The best images were selected by calculating the ratio be- 
tween peak flux and total flux of the observed star, referred to a 
circular aperture centered on the brightest pixel. This number is 
linearly related to the Strehl number - at least in the high-quality 
images we are interested in - and a quite reliable measure to find 
the frames least affected by atmospheric turbulence. 
The final image was constructed by combining the best 2% of 
all images using the Drizzle algorithm (Fruchter & Hook 2002), 
corresponding to a total integration time of 4.5 sec. The images 
were two times oversampled during drizzling, leading to a final 
pixel scale of 23.3 mas/px. The total time needed for data acqui- 



sition amounts to about 5 min per filter, including instrumental 
overheads, and approximately the same time is needed for data 
reduction. Our own version of the drizzle algorithm makes use 
of the fact that we do not expect any changes in image rotation 
or field distortion between the individual images, and thus only 
considers translations. This allows a fast and simple implemen- 
tation, and a near real-time data reduction. 

The pixel scale and image orientation were determined with 
a set of seven 'standard binary' observations, i.e., observations of 
known visual binaries (chosen from the Sixth Catalog of Orbits 
of Visual Binary Stars (Hartkopf & Mason 2006 ), orbit grades 

1 and 2) with angular separations between 076 and 2'.' 8. These 
binaries were observed several times during the night (typically 

2 standards/hour) in the same manner and with the same instru- 
mental setup as the science targets, and the final reference im- 
ages were obtained using the same reduction techniques. The 
pixel scale and rotation angle determined using these reference 
images showed a high stability within 1% and 0.5°, respectively, 
during the whole night. Although utilizing double stars as cali- 
brators for other double star measurements is far from optimal, 
we think that using a set of several stars is sufficient to determine 
sound upper limits of the measurement errors. 

2.2. HST/NICMOS 

Observations of HD 160934 with HST/NICMOS were obtained 
on June 30, 1998 (GO 7226, PI E. Becklin). Since the main aim 
of the programme was to search for faint, substellar companions, 
the coronagraphic mask in NIC2 was used. Because of this, the 
close companion to HD 160934 was not detected in the science 
data (see Lowrance et al. 2005). It is, however, detectable in the 
two acquisition frames. 

The acquisition frames were obtained with the F165M filter 
and an integration time of 0.626 s. For the analysis we assumed 
a pixel scale of 75. 10 mas/pixel. The FWHM of the PSF was 
around 2 pixel, corresponding to 150 mas, which is close to the 
diffraction limit of HST at the observing wavelength. For the 
second acquisition frame, taken at the beginning of a subsequent 
orbit, the HST guide star acquisition partially failed, resulting in 
a slightly trailed PSF. This, however, did not strongly bias the 
analysis of the data. 

Figure Q] shows a comparison of the AstraLux RG 780 and 
RG830, and the NICMOS F165M images of the HD 160934 
binary, while Table [TJ gives the dates and details of the direct 
imaging observations. 

The analysis of the binary properties is based on the pipeline 
reduced frames. Tiny Tim Version 6.3 (Krist & Hook 2004) was 
used to compute the theoretical PSF. In order to estimate the 
effect of HST "breathing" (i.e. focus changes induced by thermal 
expansion or shrinking of the optical train of HST), also slightly 
defocussed PSFs were calculated and used for the binary fitting. 

2.3. Fitting of binary parameters 

For both data sets (AstraLux and NICMOS) a binary model 
(see Bouy et al. 120031 ) was fitted to the data in order to derive 
binary separation, position angle and brightness ratio. For the 
NICMOS data, only the first acquisition frame was used. The 
slightly trailed binary PSFs due to the partial guide star acqui- 
sition failure for the second HST orbit resulted in a bias in the 
determination of the brightness ratio using a non-trailed PSF. 
While for the NICMOS data a theoretical model can serve as ref- 
erence PSF for binary fitting, this is not possible for the Lucky 
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Fig. 1. AstraLux RG780, z' and NICMOS F165M images of the binary using a logarithmic intensity scale. The field of view is 
2" x 2", and North is up and East is to the left. The binary components are clearly detectable in all three bands. Note that both the 
AstraLux and NICMOS images have been oversampled by a factor of 2 compared to the original instrumental pixel scale as part of 
the drizzle process. 



Table 2. Binary properties 



Date 


separation 


PA 




brightn. ratio 


(Filter) 


["] 


[°] 






June, 30 1998 


0.155 +0.001 


275.5 + 


0.2 


0.485 + 0.006 


(F165M) 










July, 8 2006 


0.215 ±0.002 


270.9 ± 


0.3 


0.329 + 0.051 


(RG830) 











Imaging data. Since the PSF shape depends strongly on actual 
seeing conditions, it is diffcult to predict the theoretical PSF of 
a single star in the final results. Neither is there a single star 
available in our images, which could have served as an accurate 
reference PSF. We therefore decided to use a set of eight differ- 
ent reference PSFs, generated from observations of single stars 
throughout the same night, and to use the weighted average of 
separation, brightness ratio and position angle from all eight fit 
results, using the square of the residuals as weighting parame- 
ter. While the residuals vary strongly for the different reference 
PSFs, the resulting separations and position angles show only a 
small scatter. This is not quite the case for the measured bright- 
ness ratios - here the lcr error of the weighted mean is ~16%, 
considerably higher than for the HST-based F165M brightness 
ratio. The derived fitting results can be found in Table [2] 



Table 3. Unresolved photometry of HD 160934 



Source Filter/Color mag lcr 



Konigstuhl 70cm 


V 


10.192 


+ 


0.014 




U-B 


0.947 


+ 


0.008 




B-V 


1.215 


+ 


0.005 




V-R c 


0.789 


+ 


0.004 




Rc-Ic 


0.766 


+ 


0.007 




z' 


8.820 


+ 


0.009 


Weis (11991111993* 


V 


10.28 


+ 


0.020 




U-B 


0.95 


+ 


0.015 




B-V 


1.23 


+ 


0.015 




V-R 


0.78° 


+ 


0.015 




R-I 


0.63" 


+ 


0.015 


2MASS 


J 


7.618 


+ 


0.024 




H 


6.998 


+ 


0.016 




K 


6.812 


+ 


0.020 


Hipparcos 


V 


10.29 








B-V 


1.591 


+ 


0.400 




V-I c 


2.58 


+ 


0.91 



:l The R and I-band photometry of Weis is given in the Kron system. 
Using the cubic tranformations given by Bessell and Weis t' 19871 1. 
the corresponding colors in the Cousins system are V-R f =0.78 and 
R c -I c =0.79. 



2.4. Unresolved photometry 

Available photometry of the unresolved binary covers the wave- 
length range from U- to K-band. Weis 019911 11993) has pub- 
lished Johnson-Kron UBVRI photometry. JHK magnitudes are 
contained in the 2MASS point source catalogue (Skrutskie et 
al. |20061 >. In addition, we obtained own UBV(RI)c and z' pho- 
tometry, using MPIA's Konigstuhl mountain 70cm telescope on 
September 5' /! and September 13'' 1 2006. We also performed nar- 
rowband photometry in Ha in order to derive an Ha-R color in- 
dex, but calibration of this measurement proved to be rather diffi- 
cult and inaccurate, probably due to the relatively wide passband 
of the available Ha interference filter. 

During the course of the Hipparcos mission, 96 photomet- 
ric measurements of HD 160934 were acquired. The lightcurve 
shows irregular brightness variations with semi-amplitudes of 
0.05-0.1 mag on timescales on the order of few days. Further 



evidence for the variability of HD 160934 is given by Pandey 
et al. ( 120021 1 and Henry et al. (|19951 l, suggesting amplitudes of 
«0.1mag. While Pandey et al. find a period of 43.2 days, the 
observations of Henry et al. point to a much shorter value of 
1.84 days with some uncertainties regarding the presence of a 
longer period. 

Table [3] summarises all available unresolved photometric mea- 
surements. 

2.5. Radial velocity measurements 

Radial velocity (RV) measurements of HD 160934 exist for the 
years 1995-2004. Galvez et al. (2006) published an RV curve 
based on 38 measurements, eventually leading to the classi- 
fication of HD 160934 as an SB1 spectroscopic binary with 
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an RV amplitude of Ka;7.2kms _1 . They deduce a period of 
P=6246.2318 days and an eccentricity of e=0.8028, and derive 
a spectral type of M2-M3V for the companion based on their 
mass estimates. However, the observations cover less than one 
orbit and hence only one minimum of the radial velocity. In ad- 
dition, the phase coverage is relatively sparse with only a single 
measurement for phases 0.2-0.9, and the given orbital parame- 
ters should be considered as very preliminary values. While the 
available data certainly allows to conclude that the orbit is rela- 
tively eccentric, it is this high eccentricity which makes period 
estimates without better coverage of the full RV curve unreliable. 

3. Physical properties of the HD 160934 binary 

3. 1 . Common proper motion 

The proper motion of the HD 160934 main component amounts 
to/iRA = -31.25 + 14.43mas/yr and // DEC = 59.44+ 1 1.21mas/yr. 
In the 8 years, which passed between the NICMOS and 
AstraLux observations, HD 160934 moved 250 ± 1 15 mas to the 
West, and 475 ±90 mas to the North. In the same period, the sep- 
aration between HD 160934 A and c increased by ss 60 mas, and 
the position angle decreased by ^5 deg (see Tabled. This gives 
strong evidence that both sources form indeed a physical binary. 

3.2. Photometric estimates of masses and spectral types 

Spectral types ranging from K7 (e.g., Reid et al. 1 19951 1 to M0 
(Zuckerman et al. 2004a) have been assigned to HD 160934. 
While Reid et al. ([1995 ) measured spectral indices, Zuckerman 
et al. d2004al > based their spectral typing on the optical V-K 
colour of the unresolved binary. 

We derived our own estimates of the spectral types and com- 
ponents' masses using the V-Ic color index, the V magnitude 
of the unresolved system, and the SDSS z' and F165M magni- 
tude differences of the components. We compared these values 
to the theoretical predicitions based on the models of Baraffe 
et al. (119981 abbreviated BCAH98 in the following text) for 
solar metallicity low-mass stars, searching for a mass combi- 
nation best fitting the available photometry. For our estimates 
we assumed coevality of the components. Since the published 
BCAH98 models do not directly predict z'-band magnitudes, we 
use the empirical color transforms of Jordi et al. (2006 ) to trans- 
form from RcTc to SDSS z'. Though this seems like a crude 
method and calculating appropriate z'-band magnitudes from 
model spectra should be preferred, we think that this method 
is sufficient, since we are not fitting the z' -magnitudes directly, 
but the z' -magnitude difference between the two components in- 
stead. For components of similar spectral types this approxima- 
tion is a valid way to convert between magnitude differences in 
the Johnson/Cousins and the SDSS photometric system. 

A similar approach was taken in the case of the F165M mag- 
nitude differences. Since the centers of the F165M and H pass- 
band are nearly identical, and the passbands differ only in width, 
we directly compare the F165M magnitude difference with the 
H-band magnitude differences of the models. As in the case of 
the z'-band difference, this is a sound approximation for combi- 
nations of similar spectral types. 

The BCAH98 model magnitudes were interpolated on a finer 
mass grid, and for each possible mass combination in the range 
of 0.2-1.0M o the combined V magnitude, V-Ic color index, 
and the SDSS z' and H-band brightness ratios were computed 
and compared to our measurements (in the case of the V-band 
magnitude, the comparison was made to absolute V magnitude, 



Table 4. HD160934Ac mass estimates for an assumed age of 
79 Myr 



Component 


mass [Mq] 


Teff [K] 


log L/Lq 


g 


A 


0.69 


4290 


-0.83 


4.60 


c 


0.57 


3780 


-1.23 


4.68 



based on the Hipparcos parallax). The residuals were weighted 
by the measurement errors, and the best fitting mass combina- 
tion found by determining the global minimum of the residuals. 
We performed fitting with models for different ages, namely in 
the range of 30 to 158 Myr. Minimum residuals were obtained 
with the BCAH98 model for an age of 79 Myr, and the resulting 
mass estimates are given in Table [H together with the effective 
temperatures, luminosities, and surface gravity from BCAH98. 

Using the modelled unresolved V-magnitude, we derive a 
distance module of M-m=2.81, corresponding to a distance of 
d=36.5pc or a parallax of n — 27.4mas. Compared with the di- 
rectly measured Hipparcos parallax of 40. 75+12. 06mas, this de- 
viates by 1.1 cr. Of course, since the Hipparcos parallax and its 
error were actually used in the fitting process, this photometric 
distance estimate is not an independent measurement and some- 
what circular. However, the derived values constitute a set of 
physical parameters compatible to observations within the mea- 
surement errors. 

While it is tempting to believe that our data may be sufficient 
to find not only the components' masses, but also to determine 
the age, it should be noted that especially the age is not very well 
constrained by the data. Using the BCAH98 models for 50Myr 
results in nearly equally small residuals, while the components' 
masses would then be 0.64 and 0.77 Mq, respectively. The error 
of our mass estimates should therefore assumed to be in the order 
of 0. 1 Mq. This clearly indicates that photometry-based mass es- 
timates without independent age constraints or spectroscopically 
determined effective temperatures are affected by relatively high 
uncertainties. It should also be noted that the V magnitudes of 
the BCAH98 models are known to be be rather inaccurate for 
very-low-mass stars (e.g., Allard et al. 1997), and that this may 
to some extent still be the case in the 0.7Mq regime. 

As a crosscheck, we computed the combined J, H and 
K magnitudes of the unresolved binary as predicted by the 
BCAH98 models, and compare them to the 2MASS observa- 
tions given in Table [3] The model magnitudes (using d=36.5pc) 
are J=7.59, H=6.99, K=6.86, which gives a maximum deviation 
of 2.4cr or 0.048mag in K-band. 

The derived values allow a tentative estimate of the compo- 
nents' spectral types, suggesting a combination of a K5 and an 
M0 star, which is in good agreement with the published spectral 
types of the unresolved binary (Reid et al. 1995; Zuckerman et 
al. l2004al) . 

3.3. Orbital parameters and comparison to RV data 

Between the HST and AstraLux observations, the change in pro- 
jected separation and position angle was 65mas and 4.6 degrees, 
respectively. This gives room for two possible scenarios: either 
the orbital period is considerably larger than the 8 years time dif- 
ference between the observations, or it is an integer fraction of 
it (including « 8 years as one possibility). The spectroscopically 
determined period of P~17.1yr is - at first sight - incompatible 
with this, since it would predict a difference in position angle of 
nearly 180 degrees for the direct imaging observations. A pos- 
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sible solution to this contradiction could be that HD 160934 is 
in fact a quadruple system (with the possible widely separated B 
component), and that the directly imaged companion is not iden- 
tical with the spectroscopic detection. However, as pointed out in 
Section 1231 the available RV data does not allow an unambigu- 
ous period determination. In fact, a period half the length of the 
suggested 17.1 years may still be compatible with the RV data 
and has to be confirmed or ruled out by further observations. 

If we assume a period of 8.55 years, and use our mass es- 
timates of Mi =0.69 and M2=0.57Mq, then the corresponding 
semimajor axis would be a=4.5A.U. or 0712 at a distance of 
36.5pc. Assuming an eccentricity of e=0.8, this results in a 
maximum possible separation between the two components of 
r=8.1 A.U. or 0722, respectively. This is very close to the sepa- 
ration observed with AstraLux in July 2006, which then means 
in return that the next periastron could be expected for around 
mid-2008. Further RV measurements and resolved imaging in 
the next 2-3 years are necessary to sort out the period ambiguity 
and to check whether the spectroscopically and directly detected 
companions are indeed identical. 

4. Conclusions 

By combining pre-discovery HST archive data, our own high 
angular resolution astrometry and unresolved photometry, we 
were able to derive mass and spectral type estimates for the 
HD 160934 system. These estimates are compatible with unre- 
solved 2MASS photometry, Hipparcos distance measurements, 
and existing age estimates for HD 160934 and the AB Doradus 
young moving group. We suggest that the directly imaged com- 
panion may be identical with the companion discovered by 
radial velocity measurements, but with an orbital period of 
P*8.5 years, about half the value of the published period. Further 
high angular resolution observations and radial velocity mea- 
surements in the next 2-3 years will allow to confirm or negate 
this suggestion. In the positive case, the combination of RV 
measurements and astrometry will instantly allow to compute 
a full set of orbital parameters, and to derive precise compo- 
nent masses. The knowledge of the orbit will enable the precise 
reanalysis of the Hipparcos measurements, resulting in much 
smaller errors for the parallax, distance, and distance module. 
This in return will make the HD 160934 system a valuable cali- 
brator for pre-main-sequence stellar models. 
The presented data are a good example for the power of the 
Lucky Imaging technique as a simple but effective tool for the 
discovery and follow-up of close binaries. Observations near the 
diffraction limit in the visible are possible with a fraction of 
the instrumental effort compared to adaptive optics systems or 
spaceborne observatories. 
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